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ABSTRACT 

A simple six-parameter ACDM model provides a successful fit to WMAP data. This holds both 
when the WMAP data are analyzed alone or in combination with other cosmological data. Even so, 
it is appropriate to examine the data carefully to search for hints of deviations from the now stan- 
dard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. 
The cosmological community has subjected the WMAP data to extensive and varied analyses. While 
there is widespread agreement as to the overall success of the six-parameter ACDM model, various 
"anomalies" have been reported relative to that model. In this paper we examine potential anoma- 
lies and present analyses and assessments of their significance. In most cases we find that claimed 
anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky 
simulations based on the best-fit model, one can select for low probability features of the WMAP data. 
Low probability features are expected, but it is not usually straightforward to determine whether any 
particular low probability feature is the result of the a posteriori selection or non-standard cosmology. 
Hypothesis testing could, of course, always reveal an alternative model that is statistically favored, 
but there is currently no model that is more compelling. We find that two cold spots in the map are 
statistically consistent with random cosmic microwave background (CMB) fluctuations. We also find 
that the amplitude of the quadrupole is well within the expected 95% confidence range and there- 
fore is not anomalously low. We find no significant anomaly with a lack of large angular scale CMB 
power for the best-fit ACDM model. We examine in detail the properties of the power spectrum data 
with respect to the ACDM model and find no significant anomalies. The quadrupole and octupole 
components of the CMB sky are remarkably aligned, but we find that this is not due to any single 
map feature; it results from the statistical combination of the full-sky anisotropy fluctuations. It 
may be due, in part, to chance alignments between the primary and secondary anisotropy, but this 
only shifts the coincidence from within the last scattering surface to between it and the local matter 
density distribution. While this alignment appears to be remarkable, there was no model that pre- 
dicted it, nor has there been a model that provides a compelling retrodiction. We examine claims 
of a hemispherical or dipole power asymmetry across the sky and find that the evidence for these 
claims is not statistically significant. We confirm the claim of a strong quadrupolar power asymmetry 
effect, but there is considerable evidence that the effect is not cosmological. The likely explanation 
is an insufficient handling of beam asymmetries. We conclude that there is no compelling evidence 
for deviations from the ACDM model, which is generally an acceptable statistical fit to WMAP and 
other cosmological data. 

Subject headings: cosmic background radiation - cosmological parameters - cosmology: observations 
- dark matter - early universe - instrumentation: detectors - large-scale structure 
of Universe - space vehicles - space vehicles: instruments - telescopes 
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1. INTRODUCTION 

The WMAP mission ([Bennett et all I2003aj) was de- 
signed to make precision measurements of the cosmic 
microwave background (CMB) to place constraints on 
cosmology. WMAP was specifically designed to min- 
imize systematic measurement errors so that the result- 
ing measurements would be highly reliable within well- 
determined and well-specified uncertainty levels. The 
rapidly switched and highly symmetric differential ra- 
diometer system effectively makes use of the sky as a 
stable reference load and renders most potential sys- 
tematic sources of error negligible. The spacecraft spin 
and precession paths on the sky create a highly inter- 
connected set of differential data. Multiple radiometers 
and multiple frequency bands enable checks for system- 
atic effects associated with particular radiometers and 
frequency dependencies. Multiple years of observations 
allow for checks of time-dependent systematic errors. 

The WMAP team has provided the raw time ordered 
data to the community. It has also made full-sky maps 
from these data, and these maps are the fundamental 
data product of the mission. If (and only if) the full- 
sky CMB anisotropy represented in a map is a realiza- 
tion of an isotropic Gaussian random process, then the 
power spectrum of that map contains all of the cosmolog- 
ical information. The maps and cosmological parameter 
likelihood function based on the power spectrum are the 
products most used by the scientific community. 

The WMAP team used realistic simulated time or- 
dered data to test and verify the map-making process. 
The WMAP and Cosmic Background Explorer) ( COBE) 
maps, produced by independent hardware and with sub- 
stantially different orbits and sky scannin g patterns have 
been f ound to be statistically consistent. iFreeman et ahl 
(120061) directly verifie d the fidelity of the WMAP team's 
map-makin g process (jHinshaw et al.l [20031 : iJarosik et al. 
120071 12010ft . It was indirectly verified by iWehus eTaL 
( 2009) as well. Finally, numerous CMB experiments have 
verified the WMAP sky maps over small patches of the 
sky (mainly with cross-correlation analyses), either to 
extract signal or to transfer the more precise WMAP 
calibration. 

WMAP data used alone are consistent with a six- 
parameter inflationary ACDM model that specifies the 
baryon density Qbh 2 , the cold dark matter density ft c h 2 , 
a cosmological constant 57a, a spectral index of scalar 
fluctuations n s , the optical depth to reionization r, 
and the scalar fluctuation amplitude (|Dunklev et al.l 
120091: lLarson et al] 12010ft . This ACDM model is flat, 
with a nearly (but not exactly) scale-invariant fluctua- 
tion spectrum seeded by inflation, with Gaussian ran- 
dom phases, and with statistical isotropy over the super- 
horizon sky. When WMAP data are combined with 
additional cosmological data, the ACDM model remains 
a good fit, with a narrower range of allowed parameter 
values (|Komatsu et al.l l2010'). It is remarkable that such 
diverse observations over a wide range of redshifts are 
consistent with the standard ACDM model. 

There are three major areas of future investigation: (1) 
further constrain allowed parameter ranges; (2) test the 
standard ACDM model against data to seek reliable evi- 
dence for flaws; and (3) seek the precise physical nature of 
the components of the ACDM model: cold dark matter, 



inflation, and dark energy. It is the second item that we 
address here: are there potential deviations from ACDM 
within the context of the allowed parameter ranges of the 
existing WMAP observations? 

A full-sky map T(n) may be decomposed into spherical 
harmonics Yi m as 

oo I 

T(n) = a irrXi m {n) (1) 

1=0 m=-l 

with 

a lm = J dnT(n)Y lm (n) (2) 

where n is a unit direction vector. If the CMB anisotropy 
is Gaussian-distributed with random phases, then each 
aim is independent, with a zero-mean (ai m ) — Gaussian 
distribution with 

(a>lm a l'm') = Sll'Smm'Ci (3) 

where Ci is the angular power spectrum and 8 is the 
Kronecker delta. C; is the mean variance per multipole 
moment I that would be obtained if one could take and 
average measurements from every vantage point through- 
out the universe. We have only our one sample of the 
universe, however, and its spectrum is related to the mea- 
sured aim coefficients by 

m— — / 

where ^C ; sky \ = C/ if we were able to average over an 

ensemble of vantage points. There is an intrinsic cosmic 
variance of 

In practice, instrument noise and sky masking compli- 
cate these relations. In considering potential deviations 
from the ACDM model in this paper, we examine the 
goodncss-of-fit of the Ci model to the data, the Gaus- 
sianity of the ai m derived from the map, and correlations 
between the a; m values. 

We recognize that some versions of ACDM (such as 
with multi-field inflation, for example) predict a weak de- 
viation from Gaussianity. To date, the WMAP team has 
found no such devia tions from Gaussianity. Thi s topic is 
further examined bv lKomatsu et al.l ()2010L [2009ft . Statis- 
tical isotropy is a key prediction of the simplest inflation 
theories so any evidence of a violation of rotational in- 
variance would be a significant challenge to the ACDM 
or any model based on standard inflation models. 

Anomaly claims should be tested for contamination by 
systematic errors and foreground emission, and should be 
robust to statistical methodology. Statistical analyses of 
WMAP CMB data can be complicated, and simulations 
of skies with known properties are usually a necessary 
part of the analysis. Statistical analyses must account for 
a posteriori bias, which is easier said than done. With the 
large amount of WMAP data and an enormous number 
of possible ways to combine the data, some number of low 
probability outcomes are expected. For this reason, what 
constitutes a "significant" deviation from ACDM can be 




difficult to specify. While methods to reduce foreground 
contamination (such as sky cuts, internal linear combina- 
tions, and template-based subtractions) can be powerful, 
none is perfect. Since claimed anomalies often tend to be 
at marginal levels of significance (e.g., 2 — 3<r), the resid- 
ual foreground level may be a significant consideration. 

The WMAP Science Team has searched for a num- 
ber of different potential systematic effects and placed 
quantitative upper limits on them. The WMAP team 
has extensively examined systema tic measurem e nt er - 
rors with each of its data releases: Uarosik et al.l ()2003l ) 
and IHinshaw et al.l (120031) for t he first-yea r data re- 
lease. iHinshaw et al.l (l2007fl and iPage et all (120071) for 
the three- year data release. IHinshaw et al.l (120091) for the 
five-year data release, and Uarosik et al.l (|2010t) for the 
current seven-year data release. Since those papers al- 
ready convey the extensive systematic error analysis ef- 
forts of the WMAP team, this paper focuses on the con- 
sistency of the data with the ACDM model and relies on 
the systematic error limits placed in those papers. Some 
data analysis techniques compute complicated combina- 
tions of the data where the systematic error limits must 
be fully propagated. 

This is one of a suite of papers pres enting the seven- 
year WMAP data. Uarosik et al.l (|2010D provide a discus- 
sion of the sky map s, systematic errors, and basic results. 
lLarson e~ al. (2010) derive the power spectr a and cosmo- 
logica l parameters from the WMAP data. IGold et al. I 
(2010) evaluate the foreground emission and place limits 
on the foreground contamination remain ing in the sepa- 
rated CMB data. iKomatsu et al.l (|2010f ) present a cos- 
mological interpretation of th e WMAP data comb ined 
with other cosmological data. [Wciland ct al. (201(J) an- 
alyze the WMAP observations of the outer planets and 
selected bright sources, which are useful both for an un- 
derstanding of the planets and for enabling these objects 
to serve as more effective calibration sources for CMB 
and other millimeter- wave and microwave experiments. 

This paper is organized as follows. In Section [2] we 
comment on the prominent large cold spot, nearby but 
offset from the Galactic Center region, that attracted at- 
tention when the first WMAP sky map was released 
in 2003. In Section [3] we comment on a cold spot in 
the southern sky that has attracted attention more re- 
cently. In Section [4] we assess the level of significance of 
the low value measured for the amplitude of the CMB 
quadrupole (I — 2) component. In Section [5] we discuss 
the lack of large-scale power across the sky. In Section [6] 
we assess the goodness-of-fit of the WMAP data to the 
ACDM model. In Section [7] we examine the alignment 
of the quadrupole and octupole. In Section [8] we assess 
claims of a hemispherical or dipole power asymmetry, 
and in Section [9] we assess claims of a quadrupolar power 
asymmetry. We summarize our conclusions in Section 

M 

2. COLD SPOT I, GALACTIC FOREGROUND EMISSION, 
AND THE FOUR FINGERS 

When WMAP data were first released in 2003, an im- 
age of the full sky was presented i n Galactic coordinates , 
centered on the Galactic center (|Bennett et al J [2003b) . 
Galactic emission was minimized for this image by us- 
ing an Internal Linear Combination (ILC) of WMAP 
data from independent frequency bands in such a way 



Figure 1. Top: A large coldcr-than-average region, highlighted 
by the white curve, appears prominently on the raw V-band tem- 
perature map. The full-sky map is shown with the Galactic plane 
horizontal across the center of the map with the Galactic center 
at the center of the displayed projection. Bottom: With the fore- 
ground signals strongly suppressed by the ILC technique, the high- 
lighted cold spot is seen to be at least as prominent. It is offset 
from the Galactic center in both latitude and longitude. This fact, 
combined with the fact that the clearly effective ILC foreground 
reduction does not diminish this feature, establishes that this is a 
CMB fluctuation and not a foreground effect. This feature is not 
anomalous in that simulated realizations of ACDM model skies 
routinely produce features like this. 




Figure 2. Visual inspection of the ILC map reveals four elon- 
gated valleys of cooler temperature that stretch from about the 
Galactic equator to nearly the south Galactic pole. Ridges of 
warmer-than-average temperature lie between the cooler fingers. 
These features are a consequence of large-scale power in the south- 
ern sky. It is more difficult to discern as much large-scale power 
in the northern sky. Cold Spot I is located near the northernmost 
part of one of the fingers, while Cold Spot II (within the red curve) 
is near the southernmost part of another finger. 

as to minimize signals with the frequency spectra of the 
Galactic foregrounds. The seven-year raw V-band map 
and ILC map are shown in Figure [TJ A prominent cold 
(blue) spot is seen near the center of these maps, roughly 
half of which lies within the KQ85y7 Galaxy mask. For 
the portion within the mask, the ILC process removes 
> 99% of the V-band pixel-pixel variance, while making 
almost no difference to the variance in the portion out- 
side the mask. In the ILC map, the variances in the two 
regions are nearly equal. Given its CMB-like spectrum 
and the fact that it is not centered on the Galactic cen- 
ter, this cold spot is very unlikely to be due to galactic 
foreground emission. Several years of modeling the sep- 
aration of foreground emission from the CMB emission 
continues to support the conclusion that this cold spot 
is not dominated by Galactic foreground emission, but 
rather is a fluctuation of the CMB. Further, the proba- 
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bility that a randomly located feature of this size would 
be near the Galactic center is ~ 5% and the probability 
of such a feature overlapping the Galactic plane, at any 
longitude, is much higher. This large central cold spot 
is a statistically reasonable CMB fluctuation within the 
context of the ACDM model. 

Since people are highly effective at detecting patterns, 
it is not surprising that a visual inspection of the WMAP 
sky map reveals interesting features. Four elongated cold 
(blue) fingers stretching from about the Galactic equator 
to the south Galactic pole are seen in Figure [5J There 
do not appear to be any similar fingers or features in the 
northern Galactic hemisphere aside from the northern- 
most extensions of the mostly southern fingers. Cold 
Spot I can be seen to be the northern part of one of the 
colder fingers. 

There may be a tendency to overestimate the signif- 
icance of features like Cold Spot I or the four fingers. 
It is very hard to define quantitative statistics for such 
features due to the visual nature of their identification. 
There is also an unavoidable posterior bias when us- 
ing narrowly defined statistics targeted at particular fea- 
tures seen in our sky. In any case, visual inspection of 
simulated ACDM maps often reveals large-scale features 
such as these, without requiring any underlying statisti- 
cal fluctuation. Indeed, it is the lack of these features in 
the Northern sky that may be the more unusual situa- 
tion. 

3. COLD SPOT II 

A detection of non-Gaussianity and/or phase correla- 
tions in the WMAP ai m data would be of great inter- 
est. While a detection of non-Gaussianity could be in- 
dicative of an experimental systematic effect or of resid- 
ual foregrounds, it could also point to new cosmologi- 
cal physics. There is no single preferred test for non- 
Gaussianity. Rather, different tests probe different types 
of non-Gaussianity; therefore, it is important to subject 
t he data to a v a riety of tests. 

iVielva et al.l (|2004h used a spherical Mexican hat 
wavelet (SMHW) analysis on the first year WMAP 
data to claim a detection of a non-Gaussian signal on 
a scale of a few degrees, independent of the WMAP 
obse rving frequency. A l so ap plying the SMHW ker- 
nel, iMukheriee fc Wanel (|20Q4f ) claimed to detect non- 
Gaussianity at ~ 99% significance. The signal is a pos- 
itive kurtosis in the wavelet coefficients attributed to a 
larger than expected number of 3° to 5° cold spots in 
the so uthern Galactic hemisphere. IMukheriee fc Wand 
((2001 foun d the same result f or the ILC map. Following 
up on this, iCruz et all (|2006l ) reported that the kurto- 
sis in the wavelet distribution could be exclusively at- 
tributed to a single cold spot, which we call Cold Spot 
II, in the sky map at Galactic coordinates (I — 209°, 
b = —57°), as indicated by the red curve in Figured In 
an ana lysis of the three-year WMAP data, ICruz et all 
(2007a) reported only 1.85% of their simulations devi- 
ated from the WMAP data either in the skewness or in 
the kurtosis estimator at three different angular scales, a 
2.35<7 effect. 

In replicating the SM HW approach described above, 
iZhang fc Hutererl (|2010l ) found ~ 2.8a evidence of kur- 
tosis and ~ 2.6cr evidence for a cold spot. These values 
are in reasonable agreement with the earlier findings for 



individual statistics. Zhang and Huterer then allowed 
for a range of possibilities in disk radius, spatial filter 
shape, and the choice of non-Gaussian statistic in a "su- 
perstatistic" . They found that 23% of simulated Gaus- 
sian random skies are more unusual than the WMAP 
sky. IZhang fc Hutererl (2010) also analyzed the sky maps 
with circular disk and Gaussian filters of varying widths. 
They found no evidence for an anomalous cold spot at 
any scale when compared with random Gaussian simu- 
lations. When requiring the SMHW spatial filter shape, 
15% of simulated Gaussian random skies were more un- 
usual than the WMAP sky using the constrained super- 
statistic. 

With 1.85% - 15% of random Gaussian skies more de- 
viant than WMAP data in a wavelet analysis (depend- 
ing on the marginalization of posterior choices) poten- 
tial physical interpretations have been proposed for this 
1.45cr - 2.35a effect. In theory, cold spots in the CMB 
can be produced by the integrated Sachs-Wolfe (ISW) 
effect as CMB photons traverse cosmic voids along the 
line of sight. If Cold Spot II is due to a cosmic void, it 
would have profound implications because ACDM does 
not produce voids of sufficient magnitude to explain it. 
IMota et all (|2008j) examined void formation in models 
where dark energy was allowed to cluster and concluded 
that voids of sufficient size to explain C old S pot II were 
not readily produced. Rudnic k et al.l (|2007l ) examined 
number counts and smoothed surface brightness in the 
NRAO VLA Sky Survey (NVSS) radio source data. They 
claimed a 20% - 45% deficit in the NVSS smoothed sur- 
face brightness in the directio n of Cold Spot II . How - 
ever, this claim was refuted bv lSmith fc Hutererl (|2010[) . 
who found no evidence for a deficit, after accounting 
for systematic effects and posterior choic es made in as- 
sessin g statistical significance. Further, iGranett et al.l 
(2010) imaged several fields within Cold Spot II on the 
Canada-France-Hawaii Telescope and attained galaxy 
counts that rule out a 100 Mpc radius spherical void at 
h igh significance, fin ding no evidence for a supervoid. 

ICruz et all (|2007b| ) suggested that the cold spot could 
be the signature of a topological defect in the form of a 
cosmic texture rather than an adiaba tic fluctuation. This 
sug gestion was f urthe r discussed bv lBridges et al.l (2008) 
and lCruz et al.l (|2008f ) : they estimate that a texture with 
Gfi ~ 1.5 x 10~ 6 could produce a cold spot with the ob- 
served properties. Independently, CMB power spectra 
combined with other cosmological data can be used to 
place l imits on a statis ti cal po pulation of topological de- 
fects. lUrrestilla et al.l (|2008t ) placed a 95% confidence 
upper limit of G^i < 1.8 x 10~ 6 based on Hubble con- 
stant, nucleosynthesis, and CMB (including three-year 
WMAP ) data. Textures at this level are compatible 
with the cold spot and arc neither favored nor disfavored 
by parameter fits- F ollowing the method described in 
lUrrestilla etH] (2008), we now place a power spectrum 
based 95% confidence upper limit of G/i < 1.5 x 10~ 6 us- 
ing the seven-year WMAP data, finer scale CMB data, 
and the Hubble constant. Since the new 95% CL upper 
limit derived from the power spectrum matches the value 
needed to explain the cold spot within the simple texture 
model previously discussed in the literature, that model 
is now disfavored; however, a more sophisticated model 
might still allow a texture interpretation. 
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In conclusion, there are two possible points of view. 
One is that the cold spot is anomalous and might be 
produced by a texture or other mechanism; this cannot 
be ruled out with current data. The other view is that 
the < 2.35cr (after posterior correction) statistical evi- 
dence for a cold spot feature is not compelling, and that 
the texture explanation is disfavored. Had the anomaly 
been significant at the part per million level instead of a 
part per thousand, the posterior marginalization issues 
would be moot: the feature would be considered a real 
anomaly. If real, a void explanation would have been in 
strong conflict with the ACDM model, but a population 
of topological defects that contribute at a low level to 
the CMB power spectrum would not so much falsify the 
model as provide a small modification to it. Given the 
evidence to date, the WMAP Team is of the opinion 
that there is insufficient statistical support to conclude 
that the cold spot is a CMB anomaly relative to ACDM. 
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Figure 3. Curve is a Blackwell-Rao estimate of the relative likeli- 
hood of the quadrupole power 1(1 + l)C2/27r in /iK 2 from WMAP . 
The WMAP ILC data were smoothed to 5° and the KQ85y7 mask 
was used, both degraded to res 5. The Gibbs sampling produced 
a likelihood that has been marginalized over all other multipolcs. 
The highly non-Gaussian likelihood distribution is characteristic of 
the lowest-/ multipoles. For Z > 32 the curves become nearly Gaus- 
sian. The vertical line with the label "ACDM" is the expected 
quadrupole from the full power spectrum ACDM model best fit 
to the WMAP data. The maximum likelihood of the WMAP 
measured I = 2 quadrupole is at the vertical dotted line. These 
two values are consistent to well within the 95% confidence region. 
The WMAP quadrupole is not anomalously low. 



4. THE LOW QUADRUPOLE AMPLITUDE 

The CMB quadrupole is the largest observable struc- 
ture in our universe. The magnitude of the quadrupole 
has been known to be lower than models predict since 
it was first measur ed by COBE (|Bennett et al.l Il992t 
IHinshaw et al.llT996D . It is also the large-scale mode that 
is most prone to foreground contamination, owing to the 
disk-like structure of the Milky Way; thus estimates of 
the quadrupole require especially careful separation of 
foreground and CMB emission. The ILC method of fore- 
ground suppression is especially appropriate for large- 
scale foreground removal, since the ILC's complex small- 
scale noise properties are unimportant in this context. 
With a sky cut applied, residual foreground contamina- 
tion of the quad rupole in the ILC map is determined t o 
be insignificant (jGold et al. 1 120101: Uarosik et al.ll2010D . 
In this section we reassess the statistical significance of 
the low quadrupole power in the ILC map. 



A statistical analysis of the quadrupole must account 
for the highly non-Gaussian posterior distribution of the 
low-1 (I < 32) multipoles. In this paper we use Gibbs 
sampling of the \ow-l power spe c trum to evaluate th e 
quadrupole (jO'Dwver et al.l [20041: IDunklev et aI|[2009T) . 
The WMAP ILC data are smoothed to 5° resolution 
(Gaussian, FWHM), degraded to HEALPix resolution 
level 5 (A^idc = 32), and masked with the KQ85y7 
mask. The Gibbs sampling produces a likelihood that 
has been marginalized over all other multipoles. A 
Blackwell-Rao estimator of the form of Equa tion 19 of 
Wandelt, Larson, & Lakshminaravanan (2004) is used 
to produce the WMAP quadrupole likelihood shown in 
Figure[3] The peak of the likelihood is at 200 /iK 2 , the 
median is at 430 /tK 2 , and the mean is at 1000 /iK 2 . The 
68% confidence range extends from 80 to 700 /iK 2 , and 
the 95% confidence range extends from 40 to 3200 /iK 2 . 

Figure HI shows the cumulative quadrupole distribu- 
tion derived from 300,000 Gibbs samples. The mean 
quadrupole predicted by the best-fit ACDM model lies 
at a cumulative probability of 0.824, well within the 95% 
confidence region allowed by the data. We conclude that 
the WMAP quadrupole measurement is not anomalously 
low. Further, while alternative models that predict a 
lower quadrupole will better match this specific part of 
the data, it is impossible to significantly favor such mod- 
els on the basis of quadrupole power alone. 
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Figure 4. Cumulative distribution function of the quadrupolcs 
from the Gibbs sampling based on 300,000 points. The vertical line 
is the predicted ACDM model quadrupole value. The cumulative 
probability is 0.824 where the vertical line crosses the cumulative 
distribution function . Since the expected quadrupole from the 
ACDM model is well within the 95% confidence range of the mea- 
sured quadrupole, accounting for noise and cosmic variance, we 
conclude that the measured quadrupole is not anomalously low. 



5. THE LACK OF LARGE-SCALE POWER 

The angular correlation function complements the 
power spectrum by measuring structure in real space 
rather than Fourier space. It measures the covariancc 
of pixel temperatures separated by a fixed angle, 

C(%) = (T^j) (6) 

where i and j are two pixels on the sky separated by an 
angle 9, and the brackets indicate an ensemble average 
over independent sky samples. Expanding the tempera- 
ture in spherical harmonics, and using the addition theo- 
rem for spherical harmonics, it is straightforward to show 



() 



that C(9ij) is related to the angular power spectrum by 

(C(%)> = ^ £(2i + 1) (Q) W^icosO) (7) 
I 

where (Ci) is the ensemble-average angular power spec- 
trum, Wi is the experimental window function, and Pi 
are the Legendre polynomials. If the CMB is statisti- 
cally isotropic, C(6ij) = C(9) depends only on the sep- 
aration of pixels i and j, but not on their individual 
directions. Since we are unable to observe an ensemble 
of skies, we must devise estimates of C(9) using the tem- 
perature measured in our sky. One approach is to assume 
that the CMB is ergodic (statistically isotropic) in which 
case C (9) can be estimated by averaging all temperature 
pairs in the sky separated by an angle 9, 

C{8) = (TiTj) U lJ=e (8) 

where the brackets indicate an average over directions 
i and j such that Zij = 9 (to within a bin). Another 
approach is to estimate the angular power spectrum C; 
and to compute C(9) using Equation (J7J)- 

The angular correlation function over the full-sky ILC 
map from Equation ([7} is shown in Figure [5j As can 
be seen, C{6) lies within the 95% confidence range of 
the best-fit ACDM model for all 9, as determined by 
Monte Carlo simulations. This supports the conclusion 
that there is no statistically significant lack of large-scale 
power on the full sky. 
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Figure 5. Angular correlation function of the full-sky WMAP 
ILC map is shown (heavy black curve). For comparison, the angu- 
lar correlation function for the best-fit ACDM model is also shown 
(thin black curve), along with the associated 68% and 95% confi- 
dence ranges, as determined by Monte Carlo simulations. The an- 
gular correlation function of the full-sky map is seen to be within 
the 95% confidence range of the best-fit ACDM model. This angu- 
lar correlation function was computed from the C; power spectrum, 
but is nearly indistinguishable from a pixel pair computation. Ei- 
ther way, there is no evidence of a lack of large-scale power. 

iSpergel et all (|2003|) applied the pixel-pair estimator to 
the first-year WMAP data and found an almost complete 
lack of correlated structure at angles > 60° for the sky, 
but that calculation was with a Galactic foreground cut. 
A foreground cut was made because of the concern that 
additional power from within the Galactic cut may arise 
from foregrounds. For regions outside the cut, it was 
appreciated that systematic errors and residual Galactic 
foregrounds are far more likely to add correlated power 



to the sky maps than to remove it. They quantified the 
lack of large-angular-scale power in terms of the statistic 

,1/2 

S l/2 = J C 2 {9)d cos 8 (9) 

and found that fewer than 0.15% of simulations had lower 
values of Sx/2- A low 6*1/2 value persists in later WMAP 
s ky maps. 

ICopi et al.l (|2007l) and lCopi etHI f|2009D claimed that 
there is evidence that the WMAP temperature fluctu- 
ations violate statistical isotropy. They directly com- 
puted the angular correlation function from pixel pairs, 
as in Equation (|8|). omitting from the sum pixel pairs 
where at least one pixel was within the foreground mask. 
The KQ85 foreground mask (at that time) removed 18% 
of the pixels from the full sky (now 22% for KQ85y7), 
while KQ75 removed 29%. Copi et al. found p-values of 
ps 0.03% for their computation of Sx/2, concluding that 
the data are quite improbable given the model. The ex- 
act p- value depended on the specific choice of CMB map 
and foreground mask. ICavon I (|2010D finds no frequency 
d ependence to the effect. 

Efst athiou et al.l (2010) find that the value of Sx/2 is 
sensitive to the metho d of computation. For example, 
lEfstathiou et al.l ([2010) computed the angular correla- 
tion function using the estimator 

C W = T 5> Z + VMwCi'Plicoae) (10) 
w 

where 

MlV = 2/TT ^ \ K lml'm'\ 2 (11) 

mm' 

and Kimi'm 1 is the coupling between modes (lm) and 
{I'm 1 ) induced by the sky cut, and Cv is the pseudo- 
power spectrum obtained by transforming the sky map 
into spherical harmonics on the cut sky. This estimator 
produced a significantly larger value for Sx/2 than the 
e stimator in Equat ion ©. 

lEfstathiou et al.l (|2010| ) also reconstructed the low-Z 
multipoles across the foreground sky cut region in a 
manner that was numerically stable, without an assump- 
tion of statistical isotropy. Their method relied on the 
fact that the low multipole WMAP data are signal- 
dominated and that the cut size is modest. They showed 
that the small reconstruction errors introduce no bias 
and they did not depend on assumptions of statistical 
isotropy or Gaussianity. The reconstruction error only 
introduced a small "noise" to the angular correlation 
function without changing its shape. 

The original use of a sky cut in calculating Sx/2. was 
motivated by concern for residual foregrounds in the ILC 
map. We now recognize that this precaution was un- 
necessary as the ILC foreground residuals are relatively 
small. Values of Sx/2 are smaller on the cut sky than on 
the full sky, but since the full sky contains the superior 
sample of the universe and the cut sky estimates suf- 
fer from a loss of information, cut sky estimates must 
be considered sub- opti mal. It now appears th at the 
ISpergel et all d200l and lCopi et all (I2007L [2009h low p- 
values result from both the a posteriori definition of Sx/2 
and a chance alignment of the Galactic plane with the 
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CMB signal. The alignment involves Cold Spot I and the 
northern tips of the other fingers, and can also be seen 
i n the maps that will be discussed in Section [71 

lEfstathiou et all (|2010t l corrected the full-sky WMAP 
ILC map for the es timated ISW signal from redshift z < 
0.3 as estimated bv lFrancis fc Peacock! (|2010l ) . The result 
was a substantial increase in the Si/2- Yet there is no 
large-scale cosmological significance to the orientation of 
the sky cut or the orientation of the local distribution of 
matter with respect to us; thus the result from Spergel 
et al. and Copi et al. must be influenced by a chance 
alignment of the ISW effect and a posterior statistical 
bias in the choice of statistic. 

More generally. lHaiian et~ al. (2005) applied their bipo- 
lar power spectrum technique and found no evidence for 
a violation of statistical isotropy at 95% CL for angular 
scales corresponding to multipole moments Z < 60. 

The low value of the S\/2 integral over the large-angle 
correlation function on the cut-sky results from a poste- 
rior choice of the statistic. Further, it is a sub-optimal 
statistic in that it is not computed over the full sky. 
There is evidence for a chance alignment of the Galactic 
plane cut with the CMB signal, and there is also evidence 
of a chance alignment of the primary CMB fluctuations 
with secondary ISW features from the local density dis- 
tribution. The full-sky angular correlation function lies 
within the 95% confidence range. For all of these rea- 
sons, we conclude that the large-angle CMB correlation 
function is consistent with ACDM. 

6. THE GOODNESS OF FIT OF THE STANDARD MODEL 

The power spectrum of the WMAP data alone 
places strong constraints on possible cosm ological models 
i|Dunklev et aT1[2009l: ILarson et al.ll2010h . Plots of the x 2 
per degree of freedom of the temperature-temperature 
power spectrum data relative to the best-fit ACDM 
model are shown in Figure [51 In Figure [7J the cumulative 
probability of the WMAP data given the ACDM model 
is evaluated based on simulations. All of the simulated 
skies were calculated for the same input ACDM model, 
but each result was fit separately. The WMAP sky 
is statistically compatible with the model within 82.6% 
confidence, with an uncertainty of ~ 5%. 

The x 2 can be elevated because of excess scatter within 
each multipole relative to the experimental noise vari- 
ance. It could also be elevated because of an accumu- 
lation of systematic deviations of the model from the 
data across different multipoles, such as would happen 
if a parameter value were incorrect. Therefore, despite 
an acceptable overall x 2 , we examine other aspects of 
the power spectrum data relative to the model that may 
have been masked by the inclusion of all of the data into 
a single x 2 value. We examine both of these possibilities 
below. 

To explore the l-to-l' correlation properties of Ci, we 
compute: 



(Cf ata 



-a 



bestfit 



)(c; 



data 
Z+AZ 



/^bestfit 

z+A2 



2(C i bo3tfit +A r ;) 
(2«+l)/i v .« 



2(c bcstf,t +jVi + Ai) 2 



' (2/+2A/+l)/ i 2 1 



ky,l + Al 



— y /sky, Z /sky, Z+AZ 





. 1.0 - 



0.8 - 

Data without taregiounti removal 
Data with foreground removal 

0.6 ~ Simulation 



1.4 (b) 






KQ75 3 years 
KQ75 S years 
O.B h K075 7 years 



200 400 «00 900 10M 1200 
Multipole moment / 

Figure 6. (a) x 2 P er degree of freedom of the seven-year 
temperature- temperature power spectrum data relative to the best- 
fit ACDM model. Light gray points are from 50 simulations that 
used the same ACDM model with the appropriate noise and cosmic 
variance included, where the error bars are driven by the number of 
simulations. The data and simulations were run through the same 
data analysis pipeline. These simulations were used to help fit the 
effective fraction of the sky, / sky to use for the data analysis, (b) 
The x 2 per degree of freedom is compared for the three-year, five- 
year, and seven-year maps. Small differences in the fit model have 
a negligible effect on these plots. The x 2 per degree of freedom for 
I ~ 300 has been slightly growing with additional data, while other 
multipole moment ranges are more random with additional data, 
(c) The x 2 per degree of freedom had the same WMAP data 
been taken in a reverse time order. The I < 400 region appears 
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(d) Variation of x with the choice of Galactic foreground mask 
appears random for I > 400. 
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Figure 7. Black curve is the cumulative probability of the 
WMAP temperature data based on 499 simulations. All of the 
simulations were drawn from the same ACDM model, but x 2 was 
evaluated with respect to the best-fit model for each realization. Of 
these, 412 (82.6%) had a lower x 2 than the vertical line at 1224.6. 
Thus, the WMAP power spectrum is statistically compatible with 
the model. The red curve is a \ 2 distribution with 1170 degrees of 
freedom, shown for comparison. 



Note that the power spectrum of the iFinkbcinc r et al.l 
(1999) dust map in this multipole range is < 10 fiK in 
W band, i.e., more than an order of magnitude smaller. 

Figure [8] shows the results of l-to-l' correlation calcu- 
lations of Ci for different values of I — I', calculated both 
for simulations and for the WMAP data. For the most 
part the data and simulations are in good agreement. 
The most discrepant correlations in C; are for Al = 1 
near I ~ 320 and Al = 2 near I ~ 280. 

Motivated by the outlier Al = 1 correlation at / ~ 320 
seen in Figure El we further explore a possible even I 
versus odd I effect in this portion of the power spectrum. 
(Note that this is an a posteriori selection.) We define 
an even excess statistic, Eg, which compares the mean 
power at even values of £ with the mean power at odd 
values of £, within a given £-range. It is essentially a 
measure of anticorrelation between adjacent elements of 
the power spectrum, with a sign indicating the phase of 
the anticorrelation: 
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where Cg — £(£+l)C^/2n, the superscript "obs" refers to 
the observed power spectrum, and the superscript "th" 
refers to a fiducial theoretical power spectrum used for 
normalization. From this definition, it follows that Eg > 
is an even excess and Eg < is an odd excess. The 
range of £ is small enough that the variation in Cf 1 is 
also small and convenient for normalization. We choose 
(a posteriori) to bin Eg by A£ = 50. 

An apparent positive Eg in the WMAP power spec- 
trum in the range I ~ 200 — 400 is investigated quan- 
titatively using Monte Carlo simulations. Our analysis 
is limited to 33 < £ < 600, which is the part of the 
power spectrum that is flat-weighted on the sky and 
where the Monte Carlo Apod ised Spherical Transform 
EstimatoR (jHivon et al.ll2002l ) pseudo-spectrum is used. 
Because the binning is by A£ = 50, the actual I range for 
this analysis is 50 — 599. The Monte Carlo realizations 
are CMB sky map simulations incorporating appropri- 
ate WMAP instrumental noise and beam smoothing. 
Each power spectrum, whether from observed data or 
from the Monte Carlo generator, is co-added from 861 



year-by-year cross spectra in the V and W bands, with 
weighting that accounts for the noise and the beam trans- 
fer functions. 

Figure [9] shows that an even excess of significance 
~ 2.7c is found for £ = 300 — 349. If we combine the 
two adjacent bins between £ = 250 and £ — 349, the sig- 
nificance of £i in the combined bins is ~ 2. 9c, with a 
probability to exceed (PTE) of 0.26% integrated directly 
from the Monte Carlo set (Figure H0| . However, it is im- 
portant to account for the fact that this significance level 
is inflated by the posterior bias of having chosen the £ 
range to give a particularly high value. 

We attempt to minimize the posterior bias by remov- 
ing bin selection from the Monte Carlo test. Instead of 
focusing on one bin, the revised test is based on the distri- 
bution of the maximum value of the significance £g/o~(£g) 
over all bins in each Monte Carlo realization. The 11436 
Monte Carlo realizations are split into two groups: 4000 
are used to compute the normalization a {Eg) for each 
bin, and 7436 are used to compute the distribution of 
the maximum value of £g/a(£g), giving the histogram 
that is compared to the single observed value. 

The results of the de-biased test are shown in Figure 
[TT1 In addition to the avoidance of bin selection, this test 
also incorporates negative excursions of Eg, which are ex- 
cesses of power at odd I. The test shows that the visually 
striking even excess in the £ = 300 — 350 bin is actually 
of low significance, with a PTE of 5.1% (top of Figure 
[TTj) . However, the test also shows that large excursions 
in odd-£ power are less frequent in the observed power 
spectrum than in the Monte Carlos, such that 99.3% of 
the Monte Carlos have a bin with greater odd-£ signif- 
icance than the observed data (bottom of Figure ITTj) . 
Thus there appears to be a modestly significant suppres- 
sion of odd-^ power. This effect is only slightly relieved 
by accounting for the posterior selection of the enhanced 
even excess, as seen in the bottom panel of Figure [TO 

We find no evidence for a radiometer-dependence of 
the effect. We were originally suspicious that the effect 
could arise from an interaction of the foreground mask 
with large-scale power in the map, but our simulation 
results dismissed this suspicion. 

Steps and other sharp features in the power spec- 
trum P(k) tend to be smeared out in translation to Ci 
space. For example, for the non-st andard "m e ander - 
ing" cosmological inflation model of iTve fc Xul (|2010l ) 
the scalar mode responsible for inflation meanders in a 
multi-dimensional potential. This leads to a primordial 
CMB power spectrum with complicated small-amplitude 
variations with wavenumber k. Conversion to Ci has the 
effect of a significan t amount of sm oothing (see, for ex- 
ample, Figure 1.4 of lWrightJ (|2004h ). It is not likely that 
any cosmological scenario can cause the observed odd 
excess deficit. Likewise, we are aware of no experimen- 
tal effects that could cause an odd excess deficit. We 
therefore conclude that this < 3c effect is most likely a 
statistical fluke. 

7. ALIGNED QUADRUPOLE AND OCTUPOLE 

The alignment o f the quadrupole and octupole was 
first pointed o ut bvlTegmark et al.l (120031) and later elab- 
orated on b y iSchwarz et al.| ( 20041). lLand fc Magueiiol 
(|2005aD . and lLand fc Magueiiol (|2005bl). The fact of the 
alignment is not in doubt, but the significance and im- 
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Figure 8. We compute Ai = I — V correlation properties of C\ for nearby multipoles for the WMAP data (blue) and comparable 
simulations (red). For the most part the data and simulations are in good agreement. The most discrepant correlations in Ci are for Ai = 1 
near I ~ 320 and Al ~ 2 near I ~ 280. 
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Figure 10. Even excess Ei in the observed power spectrum, in 
the I range 250 — 349, compared to a histogram of the same value 
as computed from 11,436 Monte Carlo realizations. The probabil- 
ity for a random realization to exceed the observed value (PTE) is 
0.26%, interpolated directly from the Monte Carlo cumulative dis- 
tribution. A Gaussian with the same area and standard deviation 
as the histogram is overplotted, with ±lcr indications. 



Figure 9. Top: even excess Eg in the observed power spectrum, 
in bins of A£ = 50, compared to the mean and scatter from 11,436 
Monte Carlo realizations. Bottom: Eg as in the top plot, converted 
to significance units by normalizing to the Monte Carlo scatter in 
each bin. Only the I = 250 - 299 and I = 300 - 349 bins show a 
significance greater than lcr. 

plications of the alignment are discussed here. 

Do foregrounds align t he quad rupole and dipole? 
iChiang. Naselskv fc Coles I (|2007t) conclude that 
the lowest spherical harmonic modes in the ILC 
map are significant ly cont aminated by foregrounds. 
IPark. Park, fc Gott I (|2007D find that the residual 
foreground emission in a map resulting from their 
own independent foreground analysis is not statis- 
tically important to the large-scale modes of CMB 
anisotropy. The large-scale modes of their map show 
anti-correlation with the Galactic foreground emission 
in the southern hemisphere, but they are agnostic on 



whether this is due to residual Galactic emission or 
by sim ply a matter of chance. IPark. Park, fc Gott I 
((20071) also assess the WMAP Team's ILC map and 
conclude that residual foreground emission in the ILC 
map does n ot affect the e s timat ed large-scale values 
significantly. iTegmark et al.l (|2003D also performed their 
own foreground analysis and conclude that their CMB 
map is clean enough that the lowest multipoles can 
be measured without any galaxy cut at all. They also 
point out that much of the CMB power falls within the 
Galaxy cut region, "seemingly coincidentally." In other 
words, they conclude that the residual foregrounds are 
subdominant to the intrinsic CMB signal even without 
any Galaxy cut so long as a reduced foregr ound map is 
used. Ide Oliveira-Costa k, Tegmark I (|2006l) believe that 
it is more likely that the true alignment is degraded by 
foregrounds rather than created by foregrounds. 
We determine the direction of the quadrupole by rotat- 
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Figure 11. Top: histogram of max (Ci /cr(£i )) for each Monte 
Carlo realization. The maximum is taken over 11 bins of A£ = 50, 
for i = 50 — 599. The red vertical line is the observed value of 
2.65, with a probability to exceed of 5.11%. Middle: histogram of 
max(— Ei /a(£i)), which is the greatest significance found in any i 
bin having an odd excess. Although the histograms are not normal 
distributions, the observed maximum even excess is equivalent to 
a ~ 1.6<r result, and the observed maximum odd excess of 0.34 at 
the vertical red line is equivalent to a ~ 2.5a result. Bottom: the 
histogram is the same as the one above, but restricted to a pre- 
selection of the simulations where the maximum of the significance 
of the even excess is greater than or equal to the observed even 
excess. 

ing the coordinate system until |af/| 2 + \ai-(\ 2 is max- 
imized, where £ = 2 and where ai^ m are the spherical 
harmonic coefficients. This maximizes the power around 
the equator of the coordinate system. The optimization 
is done by numerically checking the value of this quan- 
tity where the z vector of the coordinate system is ro- 
tated to the center of each half-degree pixel in a res 7 
(TVgide = 128) HEALPix pixelization. The pixel where 
this value is maximized is taken to be the "direction" 
of the quadrupole. The pixel on the opposite side of the 
sphere has the same value; we arbitrarily pick one. Using 
a similar method, a "direction" is found for the octupole, 
with £ = 3. 

The probability that I — 2 and I = 3 multipoles would 
be aligned is shown in Figure [T^l The < 1° alignment 
in our sky appears to be quite improbable based upon 
random simulations of the best-fit ACDM model. 

A resolution level 7 map has 196,608 pixels of about 
0?5 diameter. The best-fit alignment axis specifies two 
pixels directly opposite each other on the sphere. The 
probability of two axes randomly aligning in the same 
pair of pixels is then 2/196,608 = 0.001%. The probabil- 
ity of getting an alignment within 0?25 of a given axis is 
0.00095%, which is close to 0.001% above. 

In an attempt to gain insight into the alignment we 
start with the ILC full-sky temperature map. We then 
produce a map of AT 2 , which is a map of anisotropy 
power. This map is constructed by smoothing the seven- 
year ILC map by 10°, removing the mean value, and 
squaring. Since the ILC map is already smoothed to 




Alignment 

Figure 12. Probability of a greater quadrupole-octupole align- 
ment is given as a function of the alignment angle, in degrees. The 
probability given here does not account for the a posteriori selec- 
tion of a multipole alignment search, nor does it take into account 
the choice of the quadrupole and the octupole moments in partic- 
ular. 

1°, the total smoothing corresponds to a Gaussian with 
FWHM = VlO 2 + 1 2 ° = 10?. 05. We then create vari- 
ous masks to probe whether the dipole-quadrupole align- 
ment can be attributed to one or perhaps two localized 
features in the map. The edges of some of the masks 
are found from contours of the AT 2 map. The contours 
were selected by eye, from a gray-scale Mollweide projec- 
tion in an image editing program, and then converted to 
HEALPix fits files. Other trial masks were chosen more 
randomly, again to probe the sensitivity of the alignment 
to different regions of the map. 

For each mask, we take the seven-year ILC map that 
has been smoothed by 10°, zero the region inside the 
mask, and take a spherical harmonic transform. From 
the a£ m coefficients, we determine the angle between the 
quadrupole and octupole. 

Figure [T2] shows the smoothed, squared temperature 
map (in color) and the effect of various masks (in gray) 
on the quadrupole-octupole alignment. Masking Cold 
Spot I eliminates any significant alignment. However, 
keeping that region but masking other regions also sig- 
nificantly reduces the quadrupole-octupole alignment. 
The posterior selection of the particular masked regions 
is irrelevant as the point is only to demonstrate that 
no single region or pair of regions solely generates the 
< 1° alignment. Rather the high degree of quadrupole- 
octupole alignment results from the statistical distribu- 
tion of anisotropy power over the whole sky. This rules 
out single-void models, a topological defect at some sky 
position, or any other such explanation. The alignment 
behaves as one would expect if it originates from chance 
random anisotropy amplitudes and phases. The align- 
ment of the I = 2 and I = 3 multipoles is intimately con- 
nected with the large-scale cool fingers and intervening 
warm regions, discussed earlier, as can be seen in Fig- 
ure H3J Although the alignment is indeed remarkable, 
current evidence is more compatible with a statistical 
combination of full-sky data than with the dominance of 
o ne or two discrete regions . 

iFrancis fc Pea cock (2010) estimated the local (z < 0.3) 
density field from the 2MASS and SuperCOSMOS galaxy 
catalogs and used that field to estimate the ISW effect 
within this volume. Large-scale features were extrap- 
olated across the Galactic plane. The effects of radial 
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Figure 13. We smooth and square the ILC map, as described in the text, to make a map of (AT) 2 to roughly visualize the anisotropy 
power in a sky map (A). Cold Spot I is the highest power region at this level of resolution, but several other regions also contribute 
substantial power at this scale. For map (A), the quadrupole and octupole are aligned to within ~ 1°. To evaluate what regions contribute 
to this alignment, we mask selected gray areas as shown in the sky maps (B) through (J). Each map is labeled by the degree of quadrupole- 
octupole alignment remaining after the gray mask is applied. Masking Cold Spot I in (B) or (H) eliminates any significant alignment. 
However, keeping those regions and masking other regions, also breaks the alignment to a significant degree. We conclude that no single 
region or pair of regions generates the alignment. Rather, the combined power contributions over a substantial fraction of the full-sky map 
cause the high degree of quadrupole-octupole alignment. Note that the chance alignment of CMB power with the Galactic cut region 
discussed in Section \E\ is apparent in map (A). 



smearing from the photometric galaxy sample were re- 
duced by taking only three thick redshift shells with 
Az = 0.1. A linear bias was used to relate galax- 
ies to density, independent of both scale and redshift 
within each of the three shells. They estimated that the 
z < 0.3 data limit conta i ns ~ 4 0% of the total ISW sig- 
nal. iFrancis fc Peacockl (|2010f ) removed their estimate 
of the ISW effect from the WMAP map. One result 
was to raise the amplitude of the quadrupole while the 
octupole amplitude was relatively unchanged. More im- 
portantly, they reported that there remains no signifi- 
cant quadrupole- octupole alignment after t he ISW re- 
moval. With the IFrancis fc Peacockl (|2010D result, the 
quadrupole-octupole alignment shifts from an early uni- 
verse property to a statistical fluke that the secondary 
anisotropy effect from the local density distribution hap- 
pens to superpose on the primordial anisotropy in such 
a way as to align the quadrupole and octupole. 

8. HEMISPHERICAL AND DIPOLE POWER ASYMMETRY 

Claims of a dipolc or hemispherical power asymmetry 
in WMAP maps have appeared in the literature since 
the release of the first-year WMAP data, with estimates 
of statistical significance ranging up to 3. 8c. We distin- 
guish between a "hemispherical" power asymmetry, in 
which the power spectrum is assumed to change discon- 
tinuously across a great circle on the sky, and a "dipole" 
power asymmetry in which the CMB is assumed modu- 




Figure 14. I = 2 quadrupole and I = 3 octupole maps are added. 
The combined map is then shown superposed on the ILC map 
from Figure[2] Note that the quadrupole and octupole components 
arrange themselves to match the cool fingers and the warm regions 
in between. The fingers and the alignment of the I = 2 and 1 = 3 
multipoles are intimately connected. 

lated by a smooth cosine function across the sky, i.e., the 
CMB is assumed to be of the form 

T(n) mo dulatcd = (1 + W • n) T(n) unmodulated- (12) 

Previous analyses of WMAP data in the literature have 
fit for either hemispherical or dipolar power asymme- 
try, and the results are qualitatively similar: asym- 
metry is found with similar direction and amplitude 
in the two cases. However, analyses that use opti- 
mal estimators (as we do in our analysi s below) have 
all studied the dipolar modu lation (e.g., iHoftuft et al.l 
I2009t lHanson fc Lewis! I2009D . Furthermore, theoret- 
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ical attempts to obtain cosmological power asymme- 
try by a ltering the statistics of the primordial fluc- 
tuations dGordon |[2007t iDonogue. Dutta fc Ross 1 120091: 



Erickcek. Kamionkowski fc Carroll1 [2008: Erickcc k et al.1 
have all found a dipolar modulation rather than 
a hemispherical modulation. Therefore, we will con- 
centrate on the dipolar modulation, defined by Equa- 
tion (p~2|) . for the sake of better comparison with both 
early universe models, and with similar analyses in the 
literature. Unambiguous evidence for power asymmetry 
would have profound implications for cosmology. 

We revisit this analysis and conclude that this claimed 
anomaly is not statistically significant, after a posteri- 
ori choices are carefully removed from the analysis. In 
looking for a power asymmetry, the most significant issue 
is removing an arbitrary choice of scale, either specified 
explicitly by a maximum multipole I, or implicitly by 
a sequence of operations such as smoothing and adding 
extra noise that define a weighting in I. 

The firs t claimed hem i spheri cal power asymmetry ap- 
peared in lEriksen et al.l (|2004D . based on the first-year 
WMAP data, where a statistic for power asymmetry was 
constructed, and its value on high-resolution WMAP 
data was compared to an ensemble of Monte Carlo sim- 
ulations in a direct frequentist approach. They quoted 
a statistical significance of 95%-99%, depending on the 
range of I selected. The details of this analysis cont ained 
many arbitrary choices. Hansen et al.l (12004, 120091) also 
used a similar methodology. In Hansen et al.l (|2009|) . the 
significance of a 2 < I < 600 hemispherical power asym- 
metry was quoted as 99.6%. 

A second class of papers used a low-resolution 
pixel likelihood form alism to study power asymmetry. 
lEriksen et al.l (|2007[ ) used this approach to search for 
a dipole power asymmetry in low-resolution three-year 
WMAP d ata and a stat i stical significance of ~ 99% was 
claimed. Ho ftuft et al.l (|2009l ) repeated the likelihood 
analysis at somewhat higher resolution and quoted a sta- 
tistical significance of 3.5c - 3.8a for different choices of 
resolution. Although the likelihood estimator contained 
fewer arbitrary choices than the preceding class of pa- 
pers, the low-resolution framework still contained an ar- 
bitrary choice of angular scale, which may be tuned (in- 
tentionally or uninte ntionally) to spuriou sly increase sta- 
tistical significance. Hoftuft et al.l (|2009f ) introduced an 
explicit cutoff multipole Z mo d and the CMB signal was as- 
sumed to be un modulated for I > l mo d and modulated for 
I < U od . Both lEriksen et all (|2007f) and Hoftuft et all 
( 2009) downgraded the data in resolution, smoothed with 
a Gaussian window, and added extra white noise. These 
processing steps implicitly defined a weighting in I where 
the power asymmetry is estimated, and introduced arbi- 
trary choices into the analysis. 

A third approach to the analysis, based on op- 
timal quadratic estim ators, recently appeared in 
Hanson fe Lewisl |2009). In this approach, the WMAP 
data were kept at full resolution and a minimum- variance 
quadratic estimator was constructed for each of the three 
vector components of the dipole modulation w t . Han- 
son and Lewis found w 97% evidence for a dipole power 
asymmetry at 2 < I < 40, and « 99.6% evidence for 
dipole power asymmetry at 2 < I < 60. However, the 
result was strongly dependent on changing the I range, 
and quickly went away for higher I. A significant shift 



was seen between the KQ75 and KQ85 masks, and be- 
tween raw and clean maps, suggesting that foreground 
contamination was not negligible. 
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Figure 15. Probability for a Monte Carlo simulation to have a 
larger dipole modulation than the co-added V+W WMAP data, 
as measured by the statistic Ki, is shown as a function of the 
maximum multipole moment i mo< j that is assumed to be modu- 
lated. This can be interpreted as the statistical significance of 
power asymmetry for a fixed value of i mo d , if one does not account 
for possible a posteriori bias when choosing i mo( j . 

Comparing t hese methods, we find that the 

Hanson fc Lewisl (|2009f ) optimal quadratic estima- 
tor has significant advantages over other analysis 
methods that have appeared in the literature. 

1. There are no arbitrary choices (such as smooth- 
ing scale) in the optimal quadratic estimator. One 
can cither look for power asymmetry in a range of 
multipolcs 2 < I < l mo d, or over the entire range 
of angular scales measured by WMAP , and the 
estimator is uniquely determined by the minimum- 
variance requirement in each case. None of the pre- 
viously considered methods had this property. 

2. There is no need to degrade the WMAP data, or 
include processing steps such as adding extra noise, 
since the optimal quadratic estimator can be effi- 
ciently computed at full W MAP resolution usi ng 
the multigrid C" 1 code from lSmith et al.l ([20071 ). 

3. Statistical significance can be assessed straightfor- 
wardly by comparing the estimator with an en- 
semble of Monte Carlo simulations. In particular, 
maximum likelihood analyses in the literature have 
assessed significance using Bayesian evidence, but 
schemes for converting the evidence integral into a 
frequentist probability are not a sufficient substi- 
tute for true Monte Carlo simulations, which di- 
rectly give the probability for a simulation to be as 
anomalous as the data. 

For these reasons, we have studied dipole power asym- 
metry using the optimal quadratic estimator. We in- 
troduce a cutoff multipole / mo d, and assume that the 
CMB is isotropic for I > Z moc i and dipole-modulated for 
2 < I < l mo d- There is an optimal quadratic estimator 
Wi for each of the three components of the (vector) mod- 
ulation w (Equation (|12p ). and an estimator k\ for the 
(scalar) amplitude of the modulation. Implementation 
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details of the estimators are presented in the Appendix, 
where we also comment on the relation with maximum 
likelihood. 

Figure [15] shows the probability that the value of the 
dipole modulation statistic k\ is larger than for the 
WMAP data, when evaluated by Monte Carlo simula- 
tion. There are choices of / mo d where the power asymme- 
try appears to have high significance. For example, when 
we chose the KQ85y7 mask and Z mo( j = 67, the probabil- 
ity for a simulation to have a value of k\ that is larger 
than the data is 0.7%. This could be interpreted as 2.5<r 
evidence for a power asymmetry, but such an interpreta- 
tion would be inflating the statistical significance since 
the choice of / mo d is an a posteriori choice. Consider an 
analogous example for the five-year Ci power spectrum. 
The power in 67=512 is high by 3.7<7, but this is not really 
a 3.7cr anomaly. Rather, it reflects the fact that there are 
a large number of I values that could have been chosen. 

Now we seek to assess the global statistical significance 
of the power asymmetry without making any a posteri- 
ori choices. Consider the probability for a Monte Carlo 
simulation to have a larger value of k\ than the WMAP 
data, as a function of i mo d' This can be interpreted as 
the statistical significance for power asymmetry in the 
range 2 < I < Z mo d, for a fixed choice of 'mod- Let rj be 
the minimum value of the probability, which we find to 
be i] = 0.012 with the KQ75y7 mask (corresponding to 
'mod = 66), or r] = 0.007 with the KQ85y7 mask (cor- 
responding to l mo d = 67). We now assess whether this 
value of rj is anomalously low. To determine this, we 
compute rj for an ensemble of Monte Carlo simulations 
where 'mod is chosen to maximize the value of rj indepen- 
dently for each simulation. 

We perform the maximization over the range 10 < 
'mod < 132. (The results depend only weakly on this 
choice of range; we have taken the upper limit of the 
range to be twice as large as the most anomalous l mo d 
in the WMAP data.) We find that the probability that 
a simulation has a value of rj that is smaller than for 
WMAP is 13% for the KQ75y7 mask, and 10% for the 
KQ85y7 mask. 

Mot ivated by the power asymmetry, lErickcek et al.1 
(2009) presented a variation of the curvaton inflation- 
ary scenario in which the curvaton has a large-amplitude 
super-horizon spatial gradient that modulates the am- 
plitude of CMB anisotropy, thereby generating a hemi- 
spheri cal power asy mmetry that could match the CMB 
data. Hirata (2009) used high-redshift quasars to place 
a limit on the gradient in the amplitude of perturbations 
that would be caused in this scenario. Their limit ruled 
out the simplest version of this curvaton spatial gradient 
scenario. Our new CMB results, presented here, largely 
remove the initial motivation for this theory. 

We conclude that there is no significant evidence for 
an anomalous dipole power asymmetry in the WMAP 
data. 

9. QUADRUPOLAR DEPENDENCE OF THE TWO-POINT 
FUNCTION 

There is another effect, related to the dipole power 
asymmetry from Section [5] in which the two-point 
function of the CMB contains a component that 
varies as a quadrupole on the sky. Motivated by 
an anisotropic model of the early universe proposed 



by lAckerman et ah] (120071) th at predicts such a signal, 
iGroeneboom fc Eriksenl (|2009[ ) used a Gibbs sampler to 
claim "tentative evidence for a preferred direction" of 
(l,b) = (110°, 10°) in the five-year WMAP map. A 
theoretical model that predicts l arge-scale quad r upola r 
anisotropy was also proposed by iGordon et al.l ([2005) . 
A similar effect could be caused by WMAP's asymmet- 
ric beams, which were not accurately represented in this 
w ork, and an algebr a ic fact or was missing in the analysis. 
In lHanson fc Lewisl ([20091 ) . the missing algebraic factor 
was corrected and the effect was verified with high sta- 
tistical significance, using an optimal cut-sky quadratic 
estimator. Optimal estimators had previously been con- 
struct ed in the all-sky case bv | Pullen fc Kamionkowski I 
(|2007ft and lDvorkin et al.l (120081) . 

Recently. IGroeneboom et al.1 ([2010h returned another 
fit, this time including polarization, the factor correc- 
tion, examinations of beam asymmetries, noise misesti- 
mation, and zodiacal dust emission. The new claim was 
9cr evidence of the preferred direction (I, b) = (96°, 30°), 
which was quite far from the original alignment direc- 
tion claimed. The new preferred direction was toward 
the ecliptic poles, strongly suggesting that this is not 
a cosmological effect. The claimed amplitude was fre- 
quency dependent, also inconsistent with a cosmological 
effect. Zodiacal d ust emission was ruled o ut as the source 
of the alignment. lHanson fc Lewisl (|2009l ) found that the 
beam asymmetry was a large enough effect to explain 
the signal, although Groeneboom et al. reported the 
opposite conclusion. The claimed statistical significance 
of the quadrupolar power asymmetry is so high that it 
seems impossible for it to be a statistical fluke or built up 
by posterior choices, even given the number of possible 
anomalies that could have been searched for. 

We have implemented the o ptimal quadratic e stima - 
tor following the approach of Hanson fc Lewisl ()2009l ) 
and confirmed that the effect exists with high statis- 
tical significance. Rather than presenting an analy- 
sis that is tied to a particular model (either cosmo- 
logical or instrumental), we have found it convenient 
to parameterize the most general quadrupolar power 
asym metry using the lang u age o f the bipolar spectrum 
from lHaiian fc Souradeepl (|2003l ). This is reviewed in 
the Appendix. The summary is that the most general 
quadrupolar anomaly can be parameterized by two quan- 
tities Af{ M and Af_ 2 i which are Z-dependent and have 
five components corresponding to the degrees of free- 
dom of a quadrupole. If statistical isotropy holds, then 
A 2 U M = Af M 2J = 0. The special case where A 2 U M 

A 2 _ 2 l ^ corresponds to an anisotropic model in which 
the local power spectrum varies across the sky (i.e., the 
quadrupolar analogue of the dipole modulation in Equa- 
tion ([12])). The special case where A 2 U M sa -2Af M 2A ^ 
corresponds to an anisotropic model in which the local 
power spectrum is isotropic, but hot and cold spots have 
preferred ellipticity where the local magnitude and orien- 
tation varies across the sky. Thus there are two indepen- 
dent "flavors" of quadrupolar anomaly; the bipolar power 
spectrum distinguishes the two and also keeps track of 
the I dependence. A proposed model for the quadrupolar 
effect in the WMAP data can be tested by computing 
the bipolar power spectrum of the model, and comparing 
with estimates of the bipolar spectrum from data. 



14 



Figure [TC] shows the components of the bipolar power 
spectrum of the WMAP data that point along the eclip- 
tic axis (i.e., Af? and Aj® 2 l in ecliptic coordinates). A 
nonzero bipolar power spectrum is seen with high sta- 
tistical significance, even in a single bin with = 50, 
confirming the existence of a quadrupolar effect. 

We implemented a number of diagnostic tests to char- 
acterize the quadrupolar effect; our findings can be sum- 
marized as follows: 

1. Only the components of the bipolar power spec- 
trum that point in the ecliptic direction (i.e., com- 
ponents with M — in ecliptic coordinates) 
contain a statistically significant signal. The com- 
ponents with M = 1 or M = 2 are consistent with 
zero within their statistical errors, even if we sum 
over all values of I to maximize signal-to-noise. 

2. The effect is larger in W-band than V-band, which 
is inconsistent with a cosmological origin. 

3. The angular dependence of the effect shows a bump 
at the scale of the first acoustic peak (I w 220), dis- 
favoring an explanation from foregrounds or noise, 
which would not be expected to show acoustic 
peaks. 

4. If we split the optimal quadratic estimator into con- 
tributions from cross correlations between differen- 
tial assemblies (DA's) in WMAP , and auto corre- 
lations in which each DA is correlated with itself, 
then we find that the amplitude of the effect is 
consistent in the two cases, disfavoring instrumen- 
tal explanations that are not highly correlated be- 
tween channels (such as striping due to 1// noise). 

5. The bipolar power spectrum of WMAP satisfies 
Ajg ss —2A^P_ 2 t , corresponding to a model in 
which the small-scale power spectrum is isotropic, 
but the shapes of hot and cold spots are not. 
(In fact, for this reason, we have used the term 
"quadrupolar effect" in this section rather than 
"quadrupolar power asymmetry", which would 
suggest that the power spectrum is modulated. We 
favor the label "effect" over "anomaly" because it 
is only an anomaly in the absence of a plausible 
cause.) 

Given the strong ecliptic alignment and that the eclip- 
tic plane was the symmetry axis of the WMAP observa- 
tions, and the non-blackbody frequency dependence, we 
conclude that this is not a ACDM anomaly. It seems very 
likely that the observed quadrupolar effect is the result of 
incomplete handling of beam asymmetries. Beam asym- 
metry generates an instrumental bipolar power spectrum 
that is consistent with all five items above, and it is diffi- 
cult to think of any other instrumental contribution that 
satisfies these properties. However, we have not yet sim- 
ulated the effects of asymmetric beams to confirm this 
explanation. A full investigation of the effect of beam 
asymmetry is underway and preliminary indications from 
our work to date are consistent with our hypothesis. 

While a detailed explanation of the quadrupolar effect 
is pending, it is important to have as much confidence as 
possible that a large anomaly in WMAP does not bias 
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Figure 16. Quadrupolar bipolar power spectra, binned with 
A? = 50, are shown separately for V-band and W-band WMAP 
data, using the KQ75y7 mask. Only the components of the bipo- 
lar power spectra that point along the ecliptic axis are shown (i.e., 
components with M = in ecliptic coordinates). A statistically 
significant quadrupolar effect is seen, even for a single frequency 
band in a single angular bin. 

the estimated power spectrum. It is reassuring (item 5 
above) that the angle-averaged power spectrum appears 
to be statistically isotropic, suggesting that the power 
spectrum is "blind" to the effect (or, less sensitive to 
beam asymmetries, assuming that is the cause). Fur- 
thermore, if beam asymmetry does turn out to explain 
th e quadrupolar e f fect, t hen the analysis in Appendix B 
of iHinshaw et al.l (|2007t ) shows independently that the 
power spectrum bias due to beam asymmetry is small. 
Nevertheless, it is important to follow up on the studies 
to date, and we plan to do so in the future. 

10. CONCLUSIONS 

In the context of this paper we take an "anomaly" 
to refer to a statistically unacceptable fit of the ACDM 
model to the C; data, a statistically significant deviation 
of the ai m from Gaussian random phases, or correlations 
between the a/ m . We are not concerned here with the 
current uncertainty range of parameter values allowed by 
the ACDM model or with whether an alternative model 
is also an acceptable fit to the data. 

Numerous claims of WMAP CMB anomalies have 
been published. We find that there are a few valuable 
principles to apply to assess the significance of suspected 
anomalies: (1) human eyes and brains are excellent at 
detecting visual patterns, but poor at assessing proba- 
bilities. Features seen in the WMAP maps, such as 
the large Cold Spot I near the Galactic center region, 
can stand out as unusual. However, the likelihood of 
such features cannot be discerned by visual inspection 
of our particular realization of the universe. (2) Monte 
Carlo simulations are an invaluable way to determine the 
expected deviations within the ACDM model. Claims of 
anomalies without Monte Carlo simulations are necessar- 
ily weak claims. (3) Some parameters are weak discrim- 
inants of cosmology because they take on a broad range 
of values for multiple realizations of the same model. (4) 
A posteriori choices can have a substantial effect on the 
estimated significance of features. For example, it is not 
unexpected to find a 2a feature when analyzing a rich 
data set in a number of different ways. However, to as- 
sess whether a particular 2a feature is interesting, one is 
often tempted to narrow in on it to isolate its behavior. 
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That process involves a posteriori choices that amplify 
the apparent significance of the feature. 

Shortly after the WMAP sky maps became avail- 
able, one of the authors (L.P.) noted that the initials 
of Stephen Hawking appear in the temperature map, as 
seen in Figure [T7J Both the "S" and "H" are beautifully 
vertical in Galactic coordinates, spaced consistently just 
above the b — line. We pose the question, what is the 
probability of this occurrence? It is certainly infinitesi- 
mal; in fact, much less likely than several claimed cos- 
mological anomalies. Yet, we do not take this anomaly 
seriously because it is silly. The Stephen Hawking ini- 
tials highlight the problem with a posteriori statistics. 
By looking at a rich data set in multiple different ways, 
unlikely events are expected. The search for statistical 
oddities must be viewed differently from tests of pre- 
determined hypotheses. The data have the power to 
support hypothesis testing rooted in ideas that are in- 
dependent of the WMAP data. We can ask which of 
two well-posed theoretical ideas is best supported by the 
data. Much of the WMAP analysis happens in a dif- 
ferent context asking, "What oddities can I find in the 
data?" 




Figure 17. "SH" initials of Stephen Hawking are shown in the 
LLC sky map. The "S" and "H" are in roughly the same font size 
and style, and both letters arc aligned neatly along a line of fixed 
Galactic latitude. A calculation would show that the probability 
of this particular occurrence is vanishingly small. Yet, there is no 
case to made for a non-standard cosmology despite this extraordi- 
narily low probability event. It is clear that the combined selection 
of looking for initials, these particular initials, and their alignment 
and location are all a posteriori choices. For a rich data set, as is 
the case with WMAP , there are a lot of data and a lot of ways 
of analyzing the data. Low probability events are guaranteed to 
occur. The a posteriori assignment of a likelihood for a particular 
event detected, especially when the detection of that event is "op- 
timized" for maximum effect by analysis choices, does not result 
in a fair unbiased assessment. This is a recurrent issue with CMB 
data analysis and is often a tricky issue and one that is difficult to 
overcome. 

For example, no one had predicted that low-1 multi- 
poles might be aligned. Rather, this followed from look- 
ing into the statistical properties of the maps. Simu- 
lations, both by the WMAP team and others, agree 
that this is a highly unusual occurrence for the standard 
ACDM cosmology. Yet, a large fraction of simulated 
skies will likely have some kind of oddity. The key is 
whether the oddity is specified in advance. 

The search for oddities in the data is essential for test- 
ing the model. The success of the model makes these 
searches even more important. A detection of any highly 
significant a posteriori feature could become a serious 
challenge for the model. The less significant features dis- 
cussed in this paper provided the motivation for consid- 
ering alternative models and developing new analysis of 



WMAP (and soon Planck) data. The oddities have trig- 
gered proposed new observations that can further test 
the models. 

It is often difficult to assess the statistical claims. It 
may well be that an oddity could be found that mo- 
tivates a new theory, which then could be tested as a 
hypothesis against ACDM. The data support these com- 
parisons. Of course, other cosmological measurements 
must also play a role in testing new hypotheses. No 
CMB anomaly reported to date has caused the scientific 
community to adopt a new standard model of cosmol- 
ogy, but claimed anomalies have been used to provoke 
thought and to search for improved theories. 

We find that Cold Spot I does not result from Galactic 
foregrounds, but rather forms the northernmost part of 
one of four cool "fingers" in the southern sky. Its am- 
plitude and extent are not unusual for ACDM. In fact, 
structures with this nature are expected. 

We find that Cold Spot II is at the southernmost end 
of a different one of the southern fingers, and it has been 
shown not to be an anomalous fluctuation. 

We find that the amplitude of the I — 2 quadrupole 
component is not anomalously low, but well within the 
95% confidence range. 

We conclude that there is no lack of large-scale CMB 
power over the full WMAP sky. The low value of the 
^-statistic integral over the large-angle correlation func- 
tion on the cut sky results from a posterior choice of a 
sub-optimal (i.e., not full sky) statistic, S1/2, a chance 
alignment of the Galactic plane cut with CMB signal, 
and a chance alignment of primary CMB fluctuation fea- 
tures with secondary ISW features from the local density 
distribution. 

We find that the quadrupole and octupole are aligned 
to a remarkable degree, but that this alignment is not due 
to a single feature in the map or even a pair of features. 
The alignment does not appear to be due to a void, for 
example. We find that the alignment is intimately associ- 
ated with the fingers of the large-scale anisotropy visible 
in the southern sky, and it results from the statistical 
combination of fluctuations over the full sky. There is 
also evidence that the alignment is due, in part, to a co- 
incidental alignment of the primary anisotropy with the 
secondary anisotropy from the local density distribution 
through the ISW effect. At the present time the remark- 
able degree of alignment appears to be no more than a 
chance occurrence, discovered a posteriori with no mo- 
tivating theory. A new compelling theory could change 
this conclusion. 

There is a portion of the power spectrum where there is 
a marginally significant lack of odd multipole power rel- 
ative to even multipole power, but overall the WMAP 
data are well fit by the ACDM model. There is no sys- 
tematic error that we are aware of that could cause the 
even power excess, nor are there any cosmological effects 
that would do so. We conclude that the even excess is 
likely a statistical fluctuation that was found a posteriori. 
No motivating theory for this phenomenon is known. 

We find that claims of hemispherical and/or dipole 
asymmetries have suffered from a posteriori choices. Af- 
ter carrying out an analysis in a manner that avoids a 
posteriori bias, we find that the evidence for a hemispher- 
ical power asymmetry is weak. 

Evidence has been reported for a significant quadrupo- 
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lar power asymmetry that does not appear to be cos- 
mological in origin, and most likely results from an in- 
complete propagation of beam asymmetries. A careful 
analysis will be a subject of future work. 

We have examined selected claims of CMB anomalies, 
but this paper is not a comprehensive review article and 
we have not attempted to address every anomaly pa- 
per in the literature. However, we can extend our re- 
sults by recognizing that various claims of anomalies are 
not necessarily independent from those that have been 
examined. For example, there are smaller-scale conse- 
quences of the fact that the value of the quadrupole 
on our sky is low (but not anomalously low) compared 
with the maximum likelihood expected value for the 
best-fit ACDM model. The amplitudes of moderate 
to small scale hot and cold spots are expected to be 
reduced statistically relative to their predicted ampli- 
tude in the best-fit ACDM model due to the reduced 
contribution from the quadrupole component. That is, 
the temperature anisotropy in a specific spot has con- 
tributions from a range of multipoles, including the 
quadrupole contribution. Since the power in the mul- 



tipoles scales as l~ 2 , the low quadrupole value in our 
sky statistically reduces the amplitude stretch of hot and 
cold spots in the map of our sky. T his should help, a t 
least in part, to explain t h e results oflHou et al.1 (l2009h 
iMonteserin et al. I d2008h. lAvaita. Weber fc Wetterich I 
(|2010f ). and lLarson fc Wandelt I (|2004f >. Another exam- 
ple where the results from this paper may relate to a 
claimed anomaly is our report of a quadrupolar power 
asymmetry. This ma y be the effect that was detected by 
iWiaux et al. I ()2006l ) through an analysis with a second 
gaussian derivative. 
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for sharing the texture power spectrum with us. We 
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l2005h. CAMB (iLewis et al.l l2000t) . and CMBFAST 
(jSeliak fc Zaldarriagal [1996) packages. We are very 
grateful to Duncan Hanson for useful discussions. 



APPENDIX 

ESTIMATORS FOR DIPOLAR AND QUADRUPOLAR ANOMALIES 

In this appendix, we present the detailed construction of the estimators used to study dipole power asymmetry in 
Section [5] and quadrupolar dependence of the two-point function in Section |H1 

Bipolar power spectrum 

If statistical isotropy is assumed, then the two-point function of the CMB is parameterized by the power spectrum 

Kraifl( 2 m 2 ) = ( _ l) mi C^i^i^2^mi,-m 2 (Al) 

The bipolar power spectrum, introduced in lHaiian fc Souradeerjl(l200l . is a formalism for analogously parameterizing 
the two-point function if the assumption of statistical isotropy is relaxed. If we decompose the two-point function 
("< 1 mii<2m 2 ) into a sum of terms which transform under rotations with total angular momentum L, then we arrive at 
the following expansion: 

(a timi a iama ) = v/(24 + l)(2£ 2 + l)£ Aj$*>/2L+T ( * q q ) ( m\ m 2 M ) (^ 2 ) 

This equation defines the bipolar power spectrum A^f ■ (Note that our normalization and sign convention differ from 

lHajian fc Souradee o (2003); this will simplify some of the equations that follow.) 
The bipolar power spectrum has the following properties: 

1. Ajffi vanishes unless — L < M < L, the triple {l\,l%,L) satisfies the triangle inequality, and (l\ +£2 + L) is even. 

2. Under rotations, A\*f transforms as a spin-L object; under the parity operation n — > (— n), it transforms as 
AfX-^(-l) L AfX- 2 

3. Symmetry: Affi = Af£. 

4. Reality: A™* = (-l) M A^~ M . 

To get some intuition for the bipolar power spectrum, we now consider a series of increasingly complicated models 
and compute the bipolar power spectrum in each case. 

Ou r fir st trivial example will be an isotropic model with power spectrum Ce- In this case comparison with Equa- 
tion (|A2j) shows that the bipolar power spectrum is given by 

A-lxli = Ce 1 Se 1 e 2 SLoSMo (A3) 

In general, the L — component A®® of the bipolar power spectrum is equal to the angle-averaged power spectrum 
Ce. (Note that the properties above imply that the only component of the bipolar power spectrum with L — is A^.) 
Components with L > will parameterize deviations from statistical isotropy. 
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Our next example (from Section[5]) is an anisotropic model which is obtained by applying a dipolar "sky" modulation 
to an isotropic CMB with power spectrum Cf. 

T(n)=(l+ J2 ^mY 1m (H) ) T(n) iso (A4) 

V M=-l / 

where Tj SO (n) is an isotropic CMB realization. To first order in the modulation w\m, a short calculation shows that 
the bipolar power spectrum is given by 

A% = C t (A5) 
aim _ A iM _ w 1M (C e -i + C e ) 

A i-i.i - A i,i-i (47r)i/2 t - Ab > 

with all other components zero. The modulation does not change the sky-averaged power spectrum A$, but the power 
spectrum in local patches near the two poles will be different. 
A very similar example is the quadrupolar "sky" modulation: 



T(n)= 1+ w 2M Y 2 M{n)\T(h) iso (A7) 

\ M=-2 ) 

with bipolar power spectrum given by: 

A% = C t (A8) 
Air = ^ (A9) 

A 2M _ A 2M _ W2M{Cl-2 + Cj) 

A t-2,t ~ A l,l-2 - (-4^1/2 ^ A1U ^ 

In the quadrupolar case, the most general anisotropic two-point function is parameterized by two ^-dependent quantities 
{Af^ g and Ajf 1 ), in contrast to the dipole case. A short calculation shows that the quadrupolar anisotropy in the 
power spectrum is proportional to 2Ai_2,i + Au, so that a model which satisfies An — 2Ag-2,£ has a roughly 
isotropic power spectrum, even though the two-point function contains a component which transforms under rotations 
as a quadrupole. 

Our final example is the anisotropic early universe model from I Ackerman et al. (20071), ui which the initial adiabatic 
curvature fluctuation £(k.) is modulated in Fourier space as follows: 



C(k) = 



1 + U '2Af^2M(k) 



M=-2 



C(k)iso (All) 



The bipolar power spectrum of this model is 



2k 2 dk 



= J^TJ2^2M J ^A, 1 (fe)\(^)P(fc) (A12) 
where Ag(k) is the angular CMB transfer function. 

Estimators: General Construction 

In this appendix we will construct estimators for the bipolar power spectrum and related quantities. We will assume 
that the bipolar power spectrum of the CMB is a linear combination of N "template" shapes A\, A^- ■ ■ An- 

N 

~ ^(^)iife (A13) 

where the coefficients U are to be estimated from data. (The choice of template shapes will be discussed shortly.) 
We assume that our template shapes satisfy A®® = 0, i.e., the templates parameterize deviations from isotropy, not 
changes in the power spectrum. 

A lengthy but straightforward calculation shows that the minimum- variance estimator ti for the template coefficients 
which is unbiased (i.e., (ti) — U) is given by 

U[a]=F i r\£ 3 [a]-M J ) (A14) 
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where the quantity £j [a] is defined by 



£M=\ E (^)i;tV(24 + 1)(24 + i)(2L+ 1) r^oo 



l\ I. 



and the matrix and vector Afj are defined by 

i^=Cov(£i[a],£,-[a]) (A16) 
A/}[o] = <f 3 -[o]> (A17) 

where the covariance and expectation value are taken over isotropic realizations of the noisy CMB a.f m . 

The estimator in Equation (|A14j) can be specialized to measure different types of statistical isotropy by making 
different choices of template shapes (A4). For example, consider the dipole modulation, parameterized as in Equa- 
tion (fT2|) by a three- vector Wi. We assume that multipoles 2 < £ < £ mo d are modulated, and multipoles i > £ mo d are 
unmodulated. The bipolar power spectrum is given by 

»=i 

where the template shapes A4 are defined by 



(A.sLM = I (V4^/3)(Q 1 + Q 3 )F* M (ei) if |4 - *a| = 1 and * < Cod fA19) 



otherwise 

When specialized to these three template shapes, Equation (|A14[) gives the optimal quadrat ic estimator W i for t he 
three components of the modulation. (A closely related estimator has also been constructed in lDvorkin et al.l ((2008).) 
We also construct an estimator K\ for the total amplitude of the dipole modulation, irrespective of its direction, by 

*i = E^ 2 ( A2 °) 

i 

We have used this estimator in Section [S] to assess statistical significance of the dipole modulation, by comparing the 
WMAP value of K\ to an ensemble of simulations. 
The quadrupole modulation can be treated analogously. For example, in Fig. [16] we have shown estimates of A|_ 2 t 

and Afl in each bin, by taking (KWbins) template shapes, corresponding to the five components of the two bipolar 
power spectr a in e ach I bin. As another example, if an optimal estimator (unbinned in t) for the primordial modulation 
in Equation (|A11() is desired, one would take five template shapes corresponding to the five components of W2M, given 
by Equation ( |A12j) . 

Relation to likelihood formalism 

We conclude this appendix by showing how the optimal estimator is related to the maximum likelihood formalism, 
for the dipole and quadrupole cases. 

First consider the dipole modulation. The likelihood function L[a|w] for the modulation Wi, given noisy CMB data 
ai m , is given by: 

L[a\w] = (const.) x Det~ 1/2 [M(u;)SA/(w;) + N] exp f-^[M(w)SM(w) + N]' 1 ^ (A21) 
where M(w) is the operator which applies the modulation to a harmonic-space map, defined by 

(M(w)a) im = f d 2 nYl m {n) (1 + w l h l ) { ^ a e , m ,Y e , m ,(h) j (A22) 

•* \l'm' ) 

We will show that the modulation w which maximizes the likelihood L[w|a] is equal to value of the optimal quadratic 
estimator Wi defined in the previous appendix, under two approximations that will be discussed further. 

First suppose that the maximum likelihood modulation is small, so that the Taylor expansion of (logL) to second 
order in na is an accurate approximation near maximum likelihood. The Taylor expansion is given by: 

logL[a|u;] w —-Hij[a)wiWj + Gi[a]wi + (const.) (A23) 



where we have defined: 

14... \n\ = - ( 

dwidw 



n,-„}= -\-J^r) (l^LM'O (.\24) 

10=0 
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= -a^C^MiSMjC^a + a^C^iM,, S}C- 1 {M j ,S}C- 1 a 
+Tr{C- 1 M i SM j ) - -Tr^" 1 ^, S}C'- 1 {M i ,5}) 
Qi[a] =a jf C~ 1 M i SC- 1 a - Tr [C^M.S] 



(A25) 



(A26) 



and Mi = 8M(w)/dwi. Second, we assume that we can approximate the a^ m -dependent quantity 7i !3 [a] by its 
expectation value: 



We would like to compare this expression to the optimal quadratic estimator lii defined in Equation (|A14[) , in the 
special case where the template shapes are given by Equation (|A19[) . In this special case, a short calculation shows 
that (Hij) = Fij, and Gj[a] = £j[a] — A/} [a], which implies (toi)ml = Wi. 

We comment briefly on the two approximations made in this appendix, namely the second-order Taylor approxima- 
tion in Equation (|A23|) and the approximation Hij [a] rs (Hij)- By the central limit theorem, both approximations 
are expected to become accurate in the limit where the number of CMB modes which contribute to the estimators is 
large. The analyses in Sections [8] and |9] have taken wide I bins (either the £ range 2 < £ < £ mo d or a series of bins with 
A£ = 50) so these approximations should be very accurate and there should be little difference in practice between a 
likelihood estimator and the optimal quadratic estimator, although we have not tested this directly. In any case, the 
optimal quadratic estimator has several important advantages over a likelihood estimator, as described in Section [5] 
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